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ABSTRACT 



The extrasolar system OGLE-2006-BLG-109L is the first multiple-p lanet sys- 



tem to be discovered by gravitational microlensing (jGaudi et al.ll2008l ); the two 
large planets that have been detected have mass ratios, semimajor axis ratios, 
and equilibrium temperatures that are similar to those of Jupiter and Saturn; 
the mass of the host star is only O.5M , and the system is more compact than 
our own Solar system. We find that in the habitable zone of the host star, the 
two detected planets resonantly excite large orbital eccentricities on a putative 
earth-mass planet, driving such a planet out of the habitable zone. We show 
that an additional inner planet of > 0.3 M e at < 0.1 AU would suppress the 
eccentricity perturbation and greatly improve the prospects for habitability of 
the system. Thus, the planetary architecture of a potentially habitable OGLE- 
2006-BLG-109L planetary system - with two "terrestrial" planets and two jovian 
planets - could bear very close resemblance to our own Solar system. 

Subject headings: planetary systems — astrobiology 

1. Introduction 

It is possible that the OGLE-2006-BLG-109L system harbors other planets, including 
earth-like planets, that are beyond the sensitivity of the microlensing observations. Here we 
use dyna mics to constrain the prospects for an earth-like habitable planet in this system. 



We adopt iKasting et al.l (]1993l )'s definition of the habitable zone (HZ) as the region in which 
an Earth-like planet with an N 2 -C02-H 2 atmosphere would be able to maintain liquid 
oceans on its surface. For an Earth-like planet, the inner boundary of the HZ is the distance 
from the star below which the planet would experience a runaway greenhouse and lose its 
oceans, and the outer boundary is the distance from the star beyond which the planet would 
condense C0 2 clouds. For OGLE-2006-BLG-109L, the HZ is ~ 0.25-0.36 AU from the star. 
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For comparison, and using the same definition, the HZ in our solar sys tem is in the 0. 9 5-1.3 7 
AU heliocentric distance range. Both these estimates are reported in iKasting et al.l (119931 ). 
in their table III. 



2. Orbital dynamics of the known planets 



First, we consider the constraints on the orbital eccentricities of the two detected massive 
planets, because these will drive considerations for the dynamics of any habitable planets. 
For the planetary parameters obtained from the observations, it is straightforward to obtain 
upper limits on the orbital eccentricities: if one orbit is nearly circular, then e\ < 0.6 or 
e 2 J$ 0.4 will avoid crossing orbits (where subscripts 1 and 2 refer to the inner planet, 
planet b, and the outer planet, planet c, respectively); for comparable orbital eccentricities, 
the upper limit is ej.,2 < 0.3. The more st ringent limits obtai n from the resonance-overlap 
criterion to avoid strong dynamical chaos (IDuncan et al.lll989l ); this also requires ei j2 < 0.3. 



We analyzed numerically the mutual gravitational perturbations of the two planets, 
assuming coplanar and nearly circular orbits as initial conditions, in a grid over the range 



of semimajor axes spanning the observational uncertainty sho wn in Fig. 



initial conditions, we used the MERCURY6 N-body integrator (jChambers 



Tj Fo r each set of 



19991 ) to calculate 



the orbital evolution of the two planets over a 1 million year duration (which is more than 
an order of magnitude longer than the secular timescales - see below). The results show that 
the semimajor axes remain very stable and the eccentricities exhibit only small amplitude 
oscillations; even in the vicinity of the 3:1 and 2:1 mean motion resonances, the maximum 
amplitude of the eccentricity perturbation was found to be ~ 0.05. 

The time evolution of small planetary eccentricities is described well by the classical 
Laplace-Lagrange linear secular perturbation theory for the eccentricity vectors, {hi,ki} = 
ej{sin zui, cosccj}, where Wi is the longitude of periapse of the i th planet; the pertur bation is 
governed by a set of coupled linear differential equations ([Murray fc Dermottlll999f ): 



dt 



N 



{hj, kj} — ^ Aji {ki, —hi} 



(1) 



i=i 



(where N = 2 for a two planet system, but the theory is generalized to any N) . The elements 
of the coefficient matrix A are proportional to the planetary masses and also are functions 
of the planetary semimajor axes. The general solution to Eq. (TT]) is: 



A? 



{hj, kj} = ^2 E j {sin(git + A), cos^t + A)} 



(2) 



i=l 
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where the gi are the eigenvalues and the corresponding eigenvectors of A. A significant 
correction is owed to the averaged effects of the nearby first order 2:1 mean m otion resonance 
(see Fig. [T]) which changes the coefficients by ~ 10% (IMalhotra et al.lll989l ); the nearby 3:1 
mean motion resonance is of second order and its effect on the secular coefficients is negligible. 
In the following analysis, we include this correction which is important in locating the secular 
resonance accurately. For the best-fit parameters, the two secular modes for the planets have 
frequencies g± = 12.6"/yr and g 2 = 60.9"/yr. 



3. Secular perturbations in the habitable zone 

Next, we analyze the gravitational perturbations of the two known planets on the orbital 
dynamics of a putative planet in the HZ. We illustrate the effects by first considering a 
massless test particle: the secular perturb ations forced on the test particle's eccentricity 



vector by the massive planets are given by ([Murray fc Dermottlll999h 



N 

Vi 



{/i f (t), **(*)} = -Y\, — {M9it + Pi), cosfot + A)}, (3) 

where are coefficients that are proportional to the planet masses and planet eccentricities, 
and also depend upon the semimajor axes of the planets and of the test particle; go is the 
so-called free precession rate of the test particle's periastron due to the average gravitational 
quadrupole moment of the system. The amplitude, ef = \Jhf + kf, is called the "forced" 
eccentricity because it is due to the forcing by the massive planets. When g « the forced 
oscillations can become very large, and the test particle is said to be in a secular resonance. 
The forced eccentricity as a function of test particle semimajor axis for the inner 1 AU of the 
OGLE-2006-BLG-109L system is shown in Fig. [2k. There are two strong secular resonances 
in this inner zone, which we label as the v\ and the v 2 resonances, associated with the g\ and 
g 2 secular modes, respectively. One of these, the v\, is located within the habitable zone of 
this system, and its width is comparable to the width of the habitable zone. Uncertainties 
in the planets' semimajor axes cause an uncertainty in the exact location of the secular 
resonance, but we find that for most of the range of planetary semimajor axes spanning the 
observational uncertainties, the secular resonance location remains within the HZ (Fig. [T]). 

Note that in Fig. [2] (and also in Fig. @] below), we adopted small eccentricities, of a 
few percent, for the known planets; the resonant eccentricity excitation of the test particle 
is proportional to the amplitude of the g\ secular mode of the planets, so larger planet 
eccentricities would lead to proportionally larger test particle excitation. For large test 
particle eccentricities, we must consider nonlinear corrections to the secular perturbation 
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theory which will limit the resonant eccentricity excitation (jMalhotral Il998l ) . We calculate 
that the maximum eccentricity excitation with the v\ secular resonance is ~ 0.55 for ~ 
0.05, and proportionally higher for higher values of e^. We also note that the linear secular 
perturbation theory is a good approximation as long as ei j2 < 0.3; because the latter is also 
the stability limit discussed in section [21 our analysis here is valid for the allowed range 
of ex,2- Our analysis is also valid for non-coplanar planetary orbits if the non-coplanarity 
is not too large, inclination < 0.3, so that the secular p erturbations of eccentricities and 
inclinations remain de-coupled ([Murray fc Dermottl 119991 ). 



The above calculation was done for a massless test particle, but its conclusions extend 
to terrestrial mass planets too. We used numerical analysis to calculate the eccentricity 
excitation of an Earth-mass planet with semimajor axis in the habitable zone; the result, 
shown by the solid points in Fig. [2h, follows closely the analytical result for the test particle. 
For such a planet, large amplitude eccentricity oscillations occur with a period ~ 0.8-4 My, 
and for most values of the semimajor axis in the HZ, the forced eccentricity will carry 
the planet out of the HZ for some portion of the planet's year during its high eccentricity 
phase (Fig. [2b). 



4. Eliminating the secular resonance from the habitable zone 

On the face of it, the analysis above would suggest that the prospects for habitability 
of OGLE-2006-BLG-109L are dim. However, we can consider ways in which the secular 
resonance might be eliminated from the habitable zone. Towards this end, we note that the 
secular resonance location is sensitive to other masses that might be present in the system 
because these would affect the secular frequencies. In particular, additional low mass inner 
planets would change go significantly while having little effect on g\ and gi- Below, we 
calculate the mass, mi P , of an additional inner planet at semimajor axis a ip inward of the 
HZ that would eliminate the secular resonance from the HZ. We also consider solutions with 
an additional outer planet. 

Additional inner planet With just the two observed planets, the free precession frequency 
go{a) is a monotonically increasing function of a. The addition of an inner planet m; p 
has the effect of changing the functional dependence of go{a) so that it has a minimum 
somewhere between the location of the inner planet at a ip and the location, a\, of planet b. 
We analytically obtain an expression for this minimum (which is a function of m; p and di P ), 
then use the condition that min{go} > gi to obtain the relationship between m ip and a- ip . 

With three perturbing planets, the free precession frequency of a HZ planet (approxi- 
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mated as a test particle) is given by (jMalhotral 1 1 9981 ) 



miaxbf) 2 (ai) 



+ 



m*a 2 wa 



+ 



3/2 



(4) 



where G is the universal const ant of gravitation, /a,cti = a/a\,at2 = a/a 2 , and 

b^) 2 (a) is a Laplace coefficient ( Murray fc Dermott 1999 ). We substitute the leading order 



approximation, b 



(i) 

3/2 



a 



3a + 0(a 3 ), in Eq. H and find 



mfn{sf } 



15 /Gma 1 / 2 ™! + (a 1 /a 2 ) 3 m 2 
14 V a? / m* 



7 



n? 



ip 



3 mi + (ai/a 2 ) 3 m 2 V a x 



"ip 



10 



(5) 



Then, using the condition min{g } > 9i to eliminate the secular resonance from the HZ, we 
obtain 



m ip > 0.0033 M a 



-V 



for a 1D <C a r 



where a res w 



j y j ^ip "Tesj (6) 

0.3 AU is the original location of the secular resonance ( unperturbed by the 
inner planet), and we used the best-fit values of m*, m 1; m 2 , «i, a 2 from I G audi et al.l (120081 ) 
to obtain the numerical coefficient. The minimum inner planet mass as a function of its 
semimajor axis is shown in Fig. [31 We see, for example, that an inner planet of ~ 0.33 M e 
at ~ 0.1 AU can eliminate the secular resonance in the HZ. A ~ 10 M e planet at ~ 0.02 AU 
(which might be called a "hot Neptune") would also be similarly effective. 

For a more complete analysis, we relax the test particle approximation for the HZ planet 
and consider the dynamics of the four planet system fully. Such a system has four secular 
frequencies g%,i — 1-4, and additional secular resonances become possible. For example, with 
a 1 M e HZ planet and an inner planet with m lv = 0.33 M e and a; p = 0.1 AU, the original 
secular resonance is displaced well interior to the HZ, but a new secular resonance arises 
near the inner edge of the HZ. We note that our perturbation analysis is not valid when the 
inner planet is too close to the HZ, i.e., a ip is in the range 0.15-0.3 AU; the perturbations 
are not small in this case, and numerical analysis indicates strong orbital instabilities and 
poor prospects for an HZ planet in this parameter regime. 

Additional outer planet An additional outer planet is also capable of suppressing the 
habitable zone secular resonance, in a slightly different way than that of the inner planet. 
An outer planet, mass m op at orbital semimajor axis a op > a 2 , has the effect of increasing 
both the frequencies go and g\, but differentially (gi increases more than g ), so that the 
secular resonance is displaced outward. The increase in g is given by 



8g 



4 m*ai p 



(7) 
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the increase in g\ is given by a similar expression but with a\ in place of a. 

To displace the secular resonance out of the HZ, we require g\ > g for all values of a 
in the HZ. This condition yields 



4 (000 - 9w)a. 



m >Z K — ^ lu/ ^ op ~ Q 44 ( Qop 

° P ~ 3 ( Gm ,)i/ 2(a f - a s/2) " U ' 44 Mffi ll AU 

where (700 and gio are the values unperturbed by the additional planet; to evaluate the 
numerical coefficient on the right hand side, we used the best-fit values for the OGLE-1006- 
BLG-109L system, and we evaluated g 00 at the outer edge of the HZ (a = 0.36 AU). The 
forced eccentricity of a 1 M e planet in such a system, with m op = 440 M e at a op = 10 AU, 
is shown in Fig. Ht>- We note that an outer planet with a semimajor axis inward of 8.5 AU 
will excite a new secular resonance inside the HZ, therefore Eq. [8] is valid for a op > 8.5 AU. 

The minimum mass estimates for the inner and outer planets, Eqs. ([HD and (181) . are 
shown in Fig. [3j Because the outer planet's minimum mass is considerably greater than a 
Jupiter mass, and is also much greater than that of the detected planets, we consider it an 
unlikely solution. 



5. Implications 

Dynamical considerations require the OGLE-2006-BLG-109L system to harbor at least 
two additional planets in order to support a potentially habitable earth-like planet. The 
solution with a sub-earth mass inner planet, an ~ earth- mass HZ planet and the two detected 
planets has an orbital structure closely resembling that of our own Solar system. Such a 
system can also support planetesimal debris belts in dynamically stable zones at 0.8-1.8 
AU and > 8.5 AU shepherded by the major planets; these would be analogous to the Solar 
system's asteroid belt and Kuiper belt, respectively. Although the solution is not unique, it 



is rea sonable from theoretical considerations for planetary system formation (e.g. iLissauer 



19931 ). and it is interesting to note that it would represent a closer resemblance to our own 



Solar system's orbital architecture than any exo-plan etary system discovered thus far. The 



1.5 kpc distance to the system (IGaudi et al.ll2008l ) is too large for current radial velocity 



and astrometric techniques to be useful to probe this system in detail, but high- sensitivity 
measurements in the infrared may potentially be able to detect the interplanetary dust 
complex generated by any debris belts in the system. 
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Planet b semimajor axis (AU) 

Fig. 1. — Semimajor axes of the two observed planets in the OGLE-2006-BLG-109L system. 
The rectangular box represents the uncertainty in the observations. The shaded area is the 
combination of planet semimajor axes that produces a secular resonance inside the habitable 
zone. Nominal locations of significant mean motion resonances are also shown. 
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Fig. 2.— Eccentricity perturbations of a particle due to the secular effects of the two 
massive planets of OGLE-2006-BLG-109L. The shaded region is the classical habitable zone, 
(a) Forced eccentricity of a massless test particle (solid line), and of an earth-mass planet 
(points), (b) Time evolution of the pericenter and apocenter distances of an earth-mass 
planet with a semimajor axis a=0.298 AU. We used the observational best-fit semimajor 
axis values of 2.3 AU and 4.6 AU and mass values of 0.71 Mj up and 0.27 Mj up for planet b 
and c, respectively, and we assumed eccentricities of e 1 = 0.03 and e 2 = 0.05. 
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Fig. 3. — The solid lines are the analytical results given in Eqs. ([6]) and (jSJ). The dashed lines 
are the result of solving by N-body numerical analysis for the minimum planet mass required 
to ensure that the forced eccentricity of a test particle is no greater than 0.1 anywhere inside 
the HZ. The shaded region is the HZ. Shown also are the semimajor axes and masses of the 
two known planets, b and c. 
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Fig. 4. — Forced eccentricity of a 1 M® planet in the OGLE-2006-BLG-109L system with 
the addition of (a) a 0.33 M® planet at 0.1 AU (shown as a large point), (b) a 440 M® 
planet at 10 AU. Dashed lines are from secular perturbation analysis; small points are from 
N-body numerical analysis. Shaded area is the habitable zone. 



